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ABSTRACT 

An  existing  program  that  determines  the  root  locus  of  n'th  order 

/ 

polynominals  has  been  modified  to  provide  plots  of  these  loci  in  the 
complex  frequency  plane  using  a  standard  line  printer.  A  methodology 
that  combines  the  computational  capabilities  of  this  root  locus  program 
with  a  variable  scale  graphical  display  of  selectable  regions  of  the 
complex  frequency  plane  is  presented.  A  listing  of  the  Fortran  IV 
source  deck  of  the  modified  program  and  two  examples  are  included. 
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I.  INTRODUCTION 

The  stability  of  linear  closed  loop  systems  Is  effectively 

determined  through  the  use  of  root  locus  methods^* For  all  but  the 

most  trivial  cases,  the  polynomial  representing  the  system’s  characteris 

tic  equation  is  high  order.  Manual  calculation  of  the  roots  of  this 

( 2 ) 

polynomial  is  laborious  and  time  consuming.  A  Fortran  IV  program 
has  been  developed  to  expedite  this  process.  The  interpretation  of  a 
tabular  display  even  when  available  is  at  bt.it  cumbersome.  The  purpose 
of  this  memo  Is  to  describe  the  use  of  a  modified  version  of  this  basic 
program  which  provides  selectable  scaled  graphical  displays  of  the 
root  loci  as  part  of  the  tabular  output. 

Main  features  of  the  modified  program  are: 

1.  Programmed  in  Fortran  IV.  No  machine  oriented  or  object 

language. 

2.  No  complex  arithmetic. 

3.  No  special  graphical  plotting  equipment  necessary. 

4.  Order  of  polynomial  may  be  no  to  100. 

5.  Number  of  variations  of  coefficient  may  be  up  to  100. 

Main  features  of  the  graphical  display  are: 

1.  Log  plot  of  third  and  fourth  quadrants  of  complex 
frequency  plane  from  0  to  10,000  radians/second. 

2.  Linear  plots  of  selected  regions  of  the  third  and  fourth 
quadrcnts  of  the  complex  frequency  plane  with  arbitrary  scales. 
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i  I .  NATURE  OF  THE  PROBLEM 


Regardless  of  the  complexity  of  a  closed  loop  system  its  transfer 
function  can  be  reduced  +o  the  equivalent  form  shown  in  Figure  I.  For 
multiple  loop  systems,  the  Grs  and  Hfs  are  readily  expressed  as  sums 
of  products  of  polynomi al s  which  are Identi f ied  with  Individual  elements 
making  up  the  complete  system. 


Figure  I.  Linear  Closed  Loop  System 


where:  R(S) 
C(S) 
B(S) 
E(S) 


System  Input 
System  Output 
System  Feedback 
System  Error 
Process  Gain 

Control ler  Gain 


The  fractions  G |  < S ) /G2 ^ S )  and  HjCSlA^IS)  are  equivalent  transfer 
functions  of  the  system  and  are  represented  by  ratios  of  polynomials 
which  upon  expansion  can  be  put  in  the  form 
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or  In  factored  fonn 


w 


n 

ir 


mfn 

ir 

1  =  1 


(S 


+  Zu> 


<S  +  Pj) 


(2) 


where  the  a's  and  b's  are  real,  the  poles  and  zeros  t-pj  and  -Z^  can 
be  either  real  or  complex  In  conjugate  pairs  and  m^  I. 

The  linear  closed  loop  system  transfer  function  Is 


C(S) 

R(S) 


KG , ( S )  H-(S> 

P  1  2 _ 

KCKpVS)  H|(S)  +  G2(S)  h2(S) 


(3) 


The  linear  closed  loop  system's  characterl stlc  equation  which 
determines  the  roots  of  the  transfer  function's  denominator  is 


KCKpG,(S)  H j  < S )  +  G2(S)  H2(S)  =  0. 

(4) 

If  we  let  KrK  =  +  K* 

C  p 

(5) 

G, (S)  H ( ( S )  =  A(S) 

G2(S)  H2(S)  =  B(S) 
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then  equation  (4)  cap  be  written  as 

K*A(S)  +  B(S)  =  0  (6) 

and  the  linear  closed  loop  system  transfer  function  is 


CCS) 

R(S) 


KG, (S)  hL(S) 

-P  1  2 

K*A(S)  +  B(S) 


The  transfer  function  for  a  unit  step  in  R(S)  becomes 


n 

CCS)  =  Kp  tt  (S  +  2^) 

u=l _ 

m+n 

S  TT  (  S  +  p  .  ) 


(7) 


(8) 


and  the  time  response  of  the  output,  c(t)  to  the  unit  step  is  given  by 

c ( t )  =  L” 1 CC(S) □  (9) 

where  L  *(  )  indicates  the  inverse  Laplace  transform 

Expressing  equation  (8)  in  partial  fraction  form  we  have 

rrH-n 


C(S) 


K  ^  K. 

+  2-,  ^ 


( 10) 


i=l  (S+p.) 


where  Kq  and  Kj  are  the  residues  at  the  respective  poles  of  C(S). 


Speci f i ca I ly 


n 

K  =  K  tt  (Z  ) 
op  ,  u 
K  u=  I 

rrH-n 

tt  (p .  ) 
i  =  l  ' 


(II) 


Kj,  =  CK  -  (S  +  Z  )] 
r  u=  I _ 

m+n 

[S  tt  ( S  +  p .  )  ] 

i  =  i  1  S=-Pj 


(12) 
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The  inverse  transform  of  (10)is  therefore 


(13) 


If  a  pole  is  complex,  its  correspond! nc  -esidue,  ,  is  also 
complex.  In  this  case  there  also  exists  a  term  of  the  form 

—  pi  t 

K  eK  where  the  bar  denotes  complex  conjugate. 


We  then  have 


where: 


<1^  ep*+  +  Kt  epl+)  =  2|K,|e°t+  cos(wtt+*  > 


( 14) 


pr  «t  + 

P(=  ot  -  J», 


♦Jl=  tan  w£ 


(15) 


Thus  if  j  poles  are  real  and  rrH-n-j  poles  are  complex  we  arrive  at  the 
complete  time  solution  for  c(t) 


c(t)  = 


K- *  *  S 


°£  + 

2|  K  I  e  cos  (wt  t  +  4 0) 


(16) 


where 

K  =  m  +  n  -  j 
2 

As  the  value  of  K*  is  varied,  the  roots  of  the  system's 
characteristic  equation  change.  Since  -  is  the  time  constant  for  a 
single  pole  and  -  is  the  damping  coefficient  associated  with  the 
frequency  w^,  we  can  see  that  the  locus  of  Pj  and  p^  in  the  complex 
plane,  as  K*  varies,  gives  an  indication  of  the  stability  of  the  system. 
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Although  the  stability  of  the  system  for  a  unit  step  function  is 
generally  adequate  the  response  of  the  system  F(t)  to  an  arbitrary  input 
R(t)  can  be  found  by  the  Duhamel  integral 

t 

F(t)  =  J  R(u)  c(t-u)  du  (17) 

o 

where  c(t)  is  of  the  form  given  in  equations  (13)  and  (13). 

A  locus  of  roots  as  K*  varies  from  K*  .  to  K*  is  cal  led  the  root 

mm  max 

locus  of  the  system.  From  (6)  it  is  seen  that  when  is  not  present 
(no  closed  loop  system),  K*=0  and  the  roots  of  the  characteristic 
equation  are  the  roots  of  B(S).  When  Kc  is  present  (closed  loop  system) 
and  the  value  of  K  approaches  infinity,  K*  also  approaches  infinity. 

The  roots  of  the  characteristic  equation  are  the  roots  of  A(S). 


III.  GRAPHICAL  METHOD 

The  root  loci  of  the  system  are  the  trajectories  of  the  roots  of 
equation  (6)  in  the  complex  plane  as  K*  varies  from  0  to  It  will  be 
recalled  that  since  the  coefficients  of  Ihe  character! sti c  equation  are 
real  the  complex  roots  appear  in  conjugate  pairs  and  the  loci  are 
symmetric  about  the  o  axis.  Hence  investigation  of  the  half  plane 
reveals  the  nature  of  the  entire  set  of  loci.  In  our  plots  of  the 
loci,  we  choose  to  display  only  the  lower  half  on  the  complex  frequency 
plane,  i.e.,  the  third  and  fourth  quadrants  which  make  up  the  left  and 
right  hand  quadrants  respectively.  We  include  two  options  (I)  a  log 
plot  and  (2)  ar  expanded  scale  linear  plot. 

In  the  log  plot  separate  plots  are  made  of  the  left  and  right 
hand  quad^nts.  Each  is  partitioned  into  dtcades  ranging  from  0.01 
rad/sec  iv  10,000  rad/sec.  Figures  3  and  4  show  the  actual  log  plots 
for  example  I,  to  he  described  later.  Figures  5  and  6  show  two  actual 
linear  expands  for  example  I,  to  be  described  later. 
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Detailed  plotting  accuracy  is  of  no  concern  on  the  log  plots. 
Selection  of  the  proper  scale  will  provide  the  desired  accuracy  on  the 
linear  expand  plots.  The  left  hand  quadrant  coordinates  are  minus 
sigma  and  j-omega.  The  right  hand  quadrant  coordinates  are  plus 

sigma  and  j-omega.  The  dimensions  of  the  coordinates  in  the 

complex  frequency  plane  are  radians/second. 

Three  different  symbols  are  shown  in  the  complex  frequency  plane 
log  plot.  The  definition  of  these  are 

*  Roots  of  B(S) 

0  Roots  of  ACS) 

•  Roots  of  K*A(S )  +  B(S) 

Dependent  on  the  proximity  of  any  of  the  above  roots  priority  is 

given  to  plotting  first  the  roots  ACS)  then  the  roots  of  BCS)  and 
finally  the  roots  of  K*A(S)  +  BCS).  For  all  real  roots  which  are 
located  on  the  plus  sigma  or  minus  sigma  axis,  the  roots  of  ACS)  appear 
on  the  sigma  axis.  Any  other  root  that  may  occur  within  the  stepping 
Increment  of  a  root  of  ACS)  i.e.,  other  roots  of  ACS),  BCS)  or 
K*ACS)  +  BCS)  will  appear  to  the  right.  The  same  priority  of  graphical 
printout  is  maintained  for  the  complex  roots. 

Rewriting  (6)  gives 

K*A(S )  +  BCS)  =  0  (6) 

Dividing  (6)  by  BCS)  and  rearranging  gives 

l~[»K»A(S)] 

BCS) 

where  -K*A(S)  =  open  loop  transfer  function  of  equivalent  system  [TFl  . 

nrsr 

The  algebraic  sign  of  the  ratio  of  the  lowest  powered  terms  of 
open  loop  transfer  function,  ACS)  will  determine  the  regions  on  the 

+  sigma  axis  where  the  roots  olS&*A(S)  are  located.  For  a  positive 

BCG ) 

ACS)  /BCS)  the 'roots  of  I  — CTFU  open  loop  =  0  are  located  on  the  zero 
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(3) 


O 


angle  locus  and  extend  from  +  infinity  on  the  positive  sigma  axis 
to  the  largest  positive  root  of  A(S)  or  B(S).  From  this  most  positive 
real  root  the  locus  of  roots  are  located  In  the  even  numbered  intervals 
between  the  roots  of  A ( S )  and  B(S)  to  minus  infinity  on  the  sigma  axis. 
For  a  A(S)/B(S)  the  roots  of  K*A(S)+B(S)  =  0  are  located  on  the  180° 
angle  locus  and  extend  from  the  largest  positive  root  of  A ( S )  or  B(S) 
on  the  Sigma  axis  leftward  to  the  adjacent  root  and  in  the  odd 
numbered  intervals  between  the  roots  of  A(S)  and  B(S)  to  minus 
infinity  on  the  sigma  axis. 

Inspection  of  the  log  plot  shown  in  Figures  3  and  4  demonstrate 

the  real  locus  for  a  positive  A(S),  while  Figures  3  and  6  demonstrate 

B  ( $ ) 

the  real  locus  for  a  negative  ^^y  When  applying  this  method  a 

6(S) 

pencilled  line  parallel  to  the  sigma  axis  to  designate  these  closed 
loop  root  trajectories  on  +he  real  axis  is  suggested. 


The  precise  location  of  these  roots  is  given  in  the  numerically 
tabulated  output.  A  detailed  discussion  of  the  complete  input  and 
output  will  be  later.  Mention  is  made  here  only  to  indicate  that  close 
continuity  exists  between  the  numerical  and  graphical  output. 

Regions  o'  interest  on  the  complex  plane  are  obvious  from  the  log 
plots.  In  order  to  study  segments  of  the  root  loci  more  closely  a 
linear  expand  plot  of  the  complex  frequency  plane  is  used.  Similar  to 
the  log  plot,  only  one  quadrant  of  the  right  and  left  hand  planes  are 
shown.  The  scaling  of  the  linear  expand  plots  is  completely  arbitrary. 
The  regions  of  interest  are  determined  from  inspection  of  the  log 
plots.  As  many  regions  as  desired  can  be  expanded.  Figures  5  and  6 
are  one  such  linear  expand  on  a  specific  region  of  the  complex 
frequency  plane. 
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IV. 


DATA  FORMATS 


A.  Polynomial  Multiplication  and  Root  Locus  Method 

1 .  General  Description. 

a.  Computes  the  numerator  polynomial  A,  where  A  Is  the  sum 
of  the  products  of  several  sets  of  polynomials.  Each  set  contains  a 
variable  number  of  polynomials  and  a  variable  number  of  sets  form  the 
sum. 

b.  Computes  the  denominator  polynomial  B,  where  B  Is  the 
sum  of  the  products  of  several  sets  of  polynomials. 

c.  Computes  the  roots  of  A. 

d.  Computes  the  roots  of  B. 

e.  Computes  the  polynomial  K  *  A  +  B,  where  K*varles  from 

K  Initial  to  K  terminate  In  Increments  of  i  K,  or  particular  values  of 
K*may  be  chosen. 

f.  Computes  the  roots  of  the  K  *  A  +  B  polynomials. 

2.  Input. 


Descrl  ptlon 

Columns 

Data 

a. 

Identification 

1 

-  80 

Identification  of  run 

b. 

Control  card 

1 

-  10 

0  for  A  K 

( Integers) 

N  for  number  of  Input 
particular  values  of  K 

1 1 

-  20 

Number  of  polynomial 
groups  to  be  added  In  A(s) 

21 

-  30 

Number  of  polynomial  groups 
to  be  added  In  B(s) 

31 

-  40 

Problem  number 

C. 

If  Ak  is  used 

1 

-  10 

K  Initial 

1 1 

-  20 

A  K 
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2.  I nput.  (Cont.) 


Description 


Columns 


Data 


cl.  If  Ak  is  used 

21-30 

C2.  If  specific  values  of  K  are  used. 


K  terminate 


I  -  10 


K 


11-20  K2 

21-30  K3 

o 

o 


o 

o 

Kki 

o  N 

(seven  values  per  card) 


d.  Numerator  group  count 
( I ntegers) 

I  -  10  Number  of  polynomials 

in  group  I  of 
numerator 

11-20  Number  of  polynomials 

in  group  2  of 
numerator 


o 

o 


o 

Etc. 

(seven  values  per  card) 


I 
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2.  I nput  (Coni . ) 


Description  Col umns  Data 

e.  Denominator  group  count 

( i  ntegers ) 

I  -  10  Number  of  polynomials 

in  group  I  of  denominator 

11-20  Number  of  polynomials 

in  group  2  of  denominator 

0 

0 

0 

Etc. 

(seven  values  per  card) 

f.  Polynomials.  Seven  values  per  card,  with  each  polynomial 
starting  on  a  new  card.  The  first  value  of  each  polynomial  Is  an 
interger,  with  the  coefficients  in  ascending  order  in  floating  point: 


1 

-  10 

Degree  +  1  of  the 

polynomi al 

1 1 

-  20 

Constant  term 

21 

-  30 

Coeff iclent  of  X 

31 

-  40 

Coefficient  of 

Etc. 

Load  in  polynomials  In  the  order  in  which  they  appear  in  the 
fraction,  with  the  numerator  polynomials  first. 

3.  Output. 

The  output  consists  of: 

a.  Value  of  K*  initial,  increment  K*,  value  of  K*  terminate. 

b.  Number  of  polynomial  groups  added  in  A(S). 

c.  Number  of  polynomial  groups  added  in  B(S). 

d.  Number  of  polynomials  in  each  group  of  A(S). 

e.  Number  of  polynomials  in  each  group  of  B(S). 

f.  Coefficients  of  polynomials  in  each  group  of  A(S)  and 

B(S) . 
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3. 


Output  (Cont. ) 


g.  The  degree,  coef f  i cients  and  roots  of  A. 

h.  The  degree,  coefficients  and  roots  of  0. 

1.  The  degree,  coefficients,  roots,  value  of  K*  for  each 
K*  A  +  B  polynomial . 

4.  Special  Considerations. 

a.  The  maximum  order  of  A  or  B  is  100. 

b.  The  maximum  number  of  specific  values  of  K  is  100. 

B.  Root  Locus  Plot  Subroutine 

I.  General  Description.  This  is  a  subroutine  designed  to  work 
with  the  basic  polynomial  multiplication  and  r'oot  locus  program.  This 
subroutine  takes  the  roots  of  K*A  +  B  calculated  by  the  polynomial 
root  locus  program  and  plots  them  on  the  complex  frequency  plane  with 
the  fol lowing  options: 

Log  plots  of  the  third  and  fourth  quadrants  of  the  complex 
frequency  plane  from  0  to  10,000  rad/sec.  Linear  expand  plots  of  the 
third  and  fourth  quadrants  of  the  complex  frequency  plane  with  select¬ 
able  scale.  The  linear  expand  plots  supplement  the  log  plot.  The 
region  covered  by  any  one  linear  expand  plot  can  be  designated  by  the 
analyst  through  the  use  of  data  cards.  The  sigma  dimensions  are 
determined  by  choosing  a  certain  point  in  either  the  third  or  fourth 
quadrant  and  picking  a  percentage  of  this  point  to  be  the  distance 
represented  on  the  sigma  axis.  (See  linear  expand  output  of  example). 
The  J-omega  dimensions  are  determined  by  choosing  a  certain  point  on 
the  J-omega  axis  and  picking  a  percentage  of  this  point  to  be  the 
distance  represented  on  the  J-omega  axis. 
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Data 


2.  I nput. 

Description 

Control  Card  Used 
to  call  the  log 
plot  subroutine 


Col umns 
I 

I  I 
21 

31  -  35 


0  (zero)  integer 

0  (zero)  integer 

0  (zero)  integer 

I0000(ten  Thousand)  integer 


Control  Cards  Used 
to  call  the  expand 
option 

Number  of  Expands 


Desl red 

1  -  10 

I  nteger 

For  each  1  1  .tear 

1  -  10 

J -omega  specific  Real 

Expand,  data  appears 

11-20 

J -omega  percent  Real 

on  one  card 

21-30 

Sigma  spec!  f  I c  Rea  1 

31-40 

Sigma  percent  Real 

Note:  There  must  be  as  many  data  cards  as  there  are  number  of 
expands  desired. 


3.  Output. 

Depending  on  the  options  specified  the  output  consists  of: 

a.  A  list  of  all  of  the  roots  plotted  on  the  log  plot. 

b.  Two  log  plots  of  the  roots  calculated  by  the 
polynomial  m j 1 1 1 p I ication  and  root  locus  program. 

c.  The  specified  number  of  linear  expand  plots  of 
selecfable  regions  in  the  complex  frequency  plane. 

V.  EXAMPLES 

Two  examples  of  the  use  of  the  techniques  described  above  will  be 
considered.  The  first  is  the  realization  of  a  typical  control  system 
and  will  be  discussed  In  some  detail.  The  second  Is  included  simply 
to  indicate  the  results  for  a  fairly  large  order  system. 
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Example  I. 


A  typical  control  system  is  shown  in  Figure  2. 


The  equivalent  form  of  Figure  2  is  shown  in  Figure  I  where 
G j ( S )  H j (S)  =  A(S)  and  G2<S)  H2(S)  =  B(S)  and  K*  =  Kp«c. 

The  characteristic  equation  for  the  system  shown  in  Figure  2  is 

-K*[2092(S2+ I .2S-I672) (S2+3.2S+25,600)  +  491 .5S2(S2+6S+400)] 

+<0.081 )(464)(S+2. 7)  <S2+3.2S  +  25,600)S 

+(0.95)  (464)  (S  +  2.7)  (S2  +  3.2S  +  25,600) 

+S(32  +  6S  +  400)  (S2  +  3.2S  +  25,600)  =  0 

Table  I  shows  the  input  card  format  for  this  characteristic  equation. 
Lines  I  through  20  apply  to  the  polynomial  multiplication  and  root 
locus  method  and  lines  21  through  28  apply  to  the  root  locus  plot 
subrout i ne. 
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Input  for  Example  I 


0 

2 

3 

331 

0.0 

0.01 

0.3 

3 

2 

4 

4 

3 

3 

25600. 

3.2 

1. 

3 

-1672. 

1.2 

1. 

1 

-2092. 

3 

0. 

0. 

-491.5 

E  s: 
—  < 

3 

400. 

6. 

1. 

t-  on 
c  o 
O  o 

1 

.081 

-I  CL 
n 

2 

0. 

464. 

Lite 

—  CO 

3 

2 

2.7 

1. 

S3 

3 

25600. 

3.2 

1. 

_j  (_ 
-  8 

1 

.95 

g  on 

z  o 

1 

464. 

—i  5 

2 

2.7 

1. 

3 

3 

25600. 

3.2 

1 . 

2 

0. 

1. 

3 

400. 

6. 

1  . 

3 

25600. 

3.2 

1. 

0 

0 

0 

10000 

6 

100. 

100. 

-50. 

100. 

LU 

z 

100. 

100. 

100. 

100. 

CO  1- 
ZD  ZD 

82 
-1  CO 
ZD 

50. 

100. 

-50. 

100. 

50. 

100. 

50. 

100. 

h-  CO 

Q  h- 

150. 

20. 

-5. 

100. 

QZ  O 
_J 
D. 

150. 

20. 

5. 

100. 
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The  tabulated  numerical  output  is  shown  in  Tables  2  through  5. 
Table  2  mirrors  the  number  of  inputs  values  of  K*,  (0=AK  increment, 

N=  number  of  particular  input  values  of  K*),  the  number  polynomial 
groups  added  in  A(S)  and  B(S),  the  problem  number,  K*  initial,  incre¬ 
ment  K*,  K*  terminate,  the  number  of  polynomial  in  each  group  of 
A(S)  and  B(S),  and  the  coefficients  of  each  of  these  polynomials 
proceeding  from  the  left  most  term  in  the  numerator  to  the  final  term 
in  the  denominator.  Each  polynomial  group  appears  on  a  separate  sheet. 
Table  3  defines  the  equivalent  open  loop  A(S)  and  B(S)  in  single 
polynomial  and  root  form.  Table  4  defines  K*  A(S)  +  B(S)  in 
polynomial  snd  root  form  for  values  of  K*  between  K*  initial  and  K* 
terminate.  Table  5  is  a  tabulation  of  K*  A(S)  *  B(S)  roots  which 
a- -3  to  be  plotted. 

Figures  3  and  4  are  the  log  plots  for  Example  I.  The  net  sign 
for  the  lowest  order  coefficients  is  positive  as  seen  in  Table  4,  hence 
the  zero  degree  locus  is  plotted.  The  roots  of  A,S)  and  B(S)  are 

A(S ) 

39.918653 

-41.162138 

-1.7304538  +  j  145.22688 
B(S) 

-1.3440477 

-1.6000000  +  j  159.99200 

-21.1 19976  +  j  20.963084 

The  locus  of  the  K*  A(S)  +  B(S)  roots  on  the  sigma  axis  as  K* 
varys  from  0  to  0.3  in  increments  of  0.010  may  be  determined  by 
inspection.  It  appears  to  the  right  of  39.918653  and  between 
-1.3440477  and  -41.162138.  The  complex  locus  has  two  branches.  One 
branch  comes  from  the  complex  pair  at  -21.119976  +  j  20.963084  and  the 
other  from  the  complex  pair  at  -1.6000000  +  j  159.99200. 

For  K*  gains  between  0  and  0.3  these  two  branches  lie  in  the  third 
and  fourth  quadrants.  Inspection  of  this  log  plot  in  conjunction  with 
the  open  loop  roots  of  A(S)  and  B(S)  given  in  Table  4  provides  insight 
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into  the  regions  of  the  complex  plane  where  an  expanded  scale  plot  of 
the  root  locus  is  required.  The  region  bounded  by  -I00.<o<  200  and 
0<w<  200  is  shown  in  Figures  5  and  6.  The  locus  between  sigma  equal  to 
-41.162138  and— 1 . 3440477  is  shown  in  detail.  The  locus  between  sigma 
equal  to  -21.119976  +  j  20.963084  and  the  positive  sigma  axis  roots 
+  39.918653  and  plus  infinity  .  s  shown  in  detail.  The  locus  in  the 
100  -  1000  rad/second  j-omega  decade  is  shown  in  a  compressed  form  in 
this  plot.  Figures  7  and  8  are  another  set  of  expanded  scale  plots  of 
the  region  bounded  by  -100  <o<l00  and  0<wl00.  This  scale  provides  more 
sensitivity  to  that  portion  of  the  locus  which  progresses  from  -21.119976 
+  20.963084  to  the  positive  sigma  axis.  The  locus  in  the  100-1000  rad/ 
sec  j-omega  decade  is  ^ot  shown.  Figures  9  and  10  are  expanded  scale 
plots  in  the  region  bounded  by  -10  <o<  10  and  20  <w<  180.  The  roots  of 
A ( b )  and  B(S)  are  in  the  third  (left  hand)  quadrant  while  the  roots  of 
K*A(S)  +  B(S)  are  in  the  fourth  (right  hand)  quadrant.  The  migration 
of  these  roots  as  K*  is  varied  is  an  indication  of  the  requirement  for 
this  expanded  scale  plot  of  the  complex  frequency  plane. 

Example  2 

Table  6  is  the  tabulated  numerical  output  for  a  more  complex 
multi  loop  system  whose  character i st i c  equation  is  of  24'th  order 
Figures  II  and  12  are  log  plots.  The  off  axis  printing  on  the  sigma 
axis  is  shown.  The  net  sign  of  the  lowest  order  coefficients  of  A(S)/ 
B(S)  is  negative,  hence  the  180°  locus  is  plotted.  Inspection  shows 
that  the  locus  is  between  +0.0040482270  and  -0.046149160  and  then 
between  each  alternate  root  or  A(S)  or  B(S)  on  the  sigma  axis. 

Inspection  of  Figures  II  and  12  indicates  that  several  regions  of  the 
complex  frequency  plane  are  of  interest.  One  such  region  is  the  region 
bounded  by  -16.0  <o<  16.0  and  25.  <w<  175.  which  is  shown  in  Figures 
13  and  14. 
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TABLE  2 

Mirror  of  Input  t  Exomple  I .  (Contd) 
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TABLE  2 

Mirror  of  Input,  Example  I  .  (Contd) 
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TABLE  2 

Mirror  of  Input,  Example  I  .  (Contd) 
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TABLE  2 

Mirror  of  Input,  Example  I.  (Contd) 
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TABLE  4 

Closed  Loop  System  Polynomials  and  Roots,  Example  I 

POLYNOMIAL  MULTIPLICATION  ANO  ROOT  LOCUS  PROBLEM  NO.  331 
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COMPLEX  FREQUENCY  PLANE.  LEFT  HAND  QUADRANT 
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Figure  5  Graphical  Display,  Linear  Expand  Plot, Left  Hand  Quadrant,  Example  I 
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Figure  6  Graph. cal  Display,  Linear  Empond  Plot,  Right  Hond  Quadronf,  E*ampie  I. 
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Figure  8.  Graphical  Display,  Linear  Expand  Plot,  Right  Hand  Quadra 


COMPLEX  FREQUENCY  PLANE,  LEFT  HAND  QUADRANT 
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LINEAR  EXPAND  PL OT < R AD/ SEC )  .00 

Figure  9  Graphical  Display,  Linear  Expand  Plot,  Left  Hand  Quadrant,  Example  I, 
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LINEAR  EXPAND  PL  01 ( RAD/ SEC  I 

Figure  !0  Graphical  Display,  Linear  Expand  Plot,  Right  Hai.d  Quadrant ,  Example  I 
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Figure  13  Graphical  Display,  Linear  Expand  Plot,  Left  Hand  Quadrand  ,  Example  2 
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VI.  CONCLUSIONS 


The  method  described  in  this  memorandum  for  computing  linear  system 
root  loci  and  plotting  their  trajectories  on  the  complex  frequency  plane 
has  the  following  merits. 

1.  Accepts  input  in  the  form  of  sums  of  products  of 
polynomials. 

2.  Produces  a  log  plot  of  the  entire  complex  plane  in  two 
graphical  displays. 

3.  Produces  selectable  scaled  linear  plots  of  regions  in  the 
complex  frequency  plane. 

4.  Relates  the  closed  loop  gain  to  the  graphical  display. 
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% 


O  CALL  KTLOCS  > 

END  ? 

C  3 

C 

c 

C  6 

SUBROUTINE  RTLOCS  ? 

DIMENSION  WLOGI 3000  I « PHIL  INI iOOO I ,MDBL IN  130001  , PHI  PIN  (  3  0CC ) ,  rf 

•  MPDBLN!  3000 )  » 

DIMENSION  SAVE  l  ( 100* 100) »  SAVE2 ( 1 00, 100) , XA ( 2000 ) ,  XH ( 2000 )  10 

DIMENSION  X( 1001* IGA ( 100) 1 1GB  I  100), C( lOOIt D( IOC), ANSI  LOO),  11 

1  SAVE  I  ICO) , ERA! 100) ,AS( 100), PSI  100) ,A( 100) ,01 100) .ROOTRI  ICC ) ,  12 

2R00TII  100)  ,  ATM  1  00  )  ,  C K  I  1  CO )  ,CK  SI  1 00  )  U 

INTEGER  COMEN  T I  20  )  14 

INTEGER  STAK1 , ST AR2 ♦ S TAR  3 , S TAK  IS 

INTEGER  ONE, ZERO, DOT, VEE, AYE  16 

INTEGER  POINT!  130,  100), DASH, BLANK, STAR, C, SPOT! 130, 50)  1) 

REAL  MDBL  I  N ,  MP D8LN  1  •» 

CATA  ONE/1HI/,ZERO/IHO/,OOT/IH./,AYE/1HA/,VEE/1MV/  19 

DATA  PLANK/1H  /, DASH/  1H-/ ,  Q/lHI/  20 

CATA  ST  ARl/ lHO/ , STAR2/IH*/ , STAK3/ 1H« /  21 

LOGICAL  PRNT  22 

EQUIVALENCE  (  SA VE 1 1  1 , 1 ) , WLCGI  l) ),( SAVE 2 ( 1 , 1 ) ,  PHIL  INI  1 ) ) ,  23 

1  I  SAVEH  l  ,31)  ,MUBLIN(  l )  )  ,ISAVE2I  l  ,  31  )  ,PHIPLN|  1  )  )  ,  24 

2  I  SAVElll,6l),MPDBLN(l)}  2*> 

EQUIVALENCE  I  SPOT  I l, 1 ) , PO  INTI l , l ) )  26 

EQUIVALENCE! STAR l, ZERO) ,<STAR3,D0T)  21 

EQUIVALENCE! JACK  I E.KJ, I  )  ,  I  JCANN, KK , J )  28 

COMMON/ ERE EK/ A, B  29 

COMMON  /BOB/  POINT  30 

COMMON  / I NF 04 /  0, BLANK, DASH  31 

COMMON  / INE05/  SAVEl,SAVE2  32 

READ!  5,  11  )  COMENT  35 

WRITEI6,  13)  COMENT  34 

13  FORMAT! ///////////, 20X.20A4)  35 

11  FORMAT!  20A4 )  36 

D0277  I OU* l , 2000  37 

XA! I OU ) =0»  0  38 

277  XBI I0U)=0.0  39 

DO  401  KK  = l , 100  40 

DO  401  KJ=l, 130  41 

401  POINTIKJ.KK )*BLANK  42 

DO  402  KK  =  1,1  30  43 

POINTIKK, 51 )=DASH  44 

402  POINTIKK, 50)=DASH  45 

DO  403  KK-1,100  46 

00  403  K J  =  34,  125,13  4  7 

403  PO I  NT ( K J , KK ) =Q  48 

DO  404  KJ=34, 124  49 

DO  404  KK  =  7 ,  93,  7  50 

IF(KK,EQ.49)  KK=5B  51 

404  POINT!KJ,KK)=DASH  52 

DO  405  KK=  7 ,  93,  7  53 

IF1KK.EQ.49)  KK  =  5B  54 

405  POINT (25, KK)= ONE  55 

DO  406  KK=26» 28  56 

406  POINTIKK,  7)=ZER0  57 

DO  407  KK=26, 27  5H 

407  POINTIKK, 14 )=ZERO  59 

POINT! 26,21 )=ZERO  oO 


61 


00  408  KK-35,65,7  61 

IFfKK.CQ.49)  KK-58  62 

408  POINT! 24, KK )*00T  63 

?0INT(2S,42I-ZER0  64 

POINT (25, 581-ZERO  65 

POINT (26*42) »ONE  66 

POINT(26,58)-ONE  67 

POINT(26,7?'-ZERO  68 

POINT(26,79)»ZERO  69 

POINT! 26, 861-ZERO  70 

POINTI27, 861-ZERO  71 

00  409  KK-26,20  72 

409  POINT! KK,93)*ZER0  73 

00  420  KK>1,100  74 

IF(KK.EQ.4)  KK-98  76 

420  POINT! 26*  KK ) -G  76 

POINT! 26f4)-VEE  77 

POINT (26, 971-AYE  70 

00  7681  JACKIE-1,100  79 

00  7681  JOANN-1,  100  bO 

SAVE!! JACKIE, J0ANN1-0.0  81 

7681  SAVE2! JACKIE, JOANNI-O.O  02 

NO-1  83 

12  00  5  I  -  1,100  04 

A!  I)  -  0.0  05 

5  B( I  I  -  0.0  86 

MI-0  07 

JZO-l  88 

Kl-1  89 

III-l  90 

REAO  (5,101)  N, I  A, IB, IPROB  91 

140  FORMAT! 10X, 7HDELTAK-, 110, 5X , 19HP0LY. ADDED  IN  A I S ) - ,  1 10 , 5X , 19HP0L Y.  92 

♦  ADDED  IN  B(S)-,I10, 5X,9HPR0B.N0. », 1 10 )  93 

IF(IB.NE.O)  60  TO  10  94 

IF!  IA.NE.0 )  GO  TO  10  95 

IF!  N.NE.O )  GO  TO  10  96 

IF! IPROB. EQ.10000)  GO  TO  1016  97 

101  FORMAT  (7110)  98 

10  WRITEI6.1)  99 

FORKl-O.O  100 

F0RK2  -  0.0  101 

F0RK3  -  0.0  102 

MRITE  (6  , 102 ) IPROB  103 

102  FORMAT  (1HI,  9X , 40HP0LYN0MI AL  MULTIPLICATION  AND  ROOT  L0CU5,44X, 11  104 

1HPR0BLEM  NO., 15//)  105 

IF  (N)  21,  20,  21  106 

20  READ( 5, 103)Y,OY,YT  107 

MRITE!6,140)N,IA, IB, IPROB  108 

WRITEI6* 141 )V,OV, YT  109 

141  FORMAT! 1H0,5X, 10HK- INITIAL-F20.10,5X, 1 2H INCREMENT  K- , F20. 10, 5X,  1 2H  110 

♦K-TERMINATE*,F20.10)  111 

103  FORMAT  I4F10.0)  112 

00  GO  TO  25  113 

CO  21  REAOI5,104)(X! I ) , I-1,N)  114 

Is-  WRITE!  6,  142)  ( I,X(  I ),  I  —  1 ,  N )  115 

<5 104  FORMAT!  7E10.0)  116 

“25  READ  ! 5  ,105)(!GA  (I),  I  «  1,  IA)  117 

NRITE( 6, 143 )  ( I , IGA( I )• I— 1, IA)  118 

143  FORMAT! 10X, 24HNUMBER  OF  POLY.  IN  GROUP,  !3,14H  OF  NUMERATOR- , I 10 )  119 

142  FORMAT! 10X,2HX(, 13, 2H)«,F20. 10)  120 


62 


READ  (5  «  1  0*>  >  (  I  GB  (  J  )  «  J=  It  IB) 

WR I T  E ( 6  » 144 )  ( J, IGB( J),  J*l, IB) 

144  FORMAT!  10X,24HNUMBER  OF  POLY.  IN  GROUP,  I3,16H  OF  OENOM  INA  I  OR  =  ,  I  lo 

*) 

105  FORMAT ( 7110) 

300  DO  15  I  =  l,  ICO 
15  SAVE! I )  =  0.0 
ISOEG  =  0 
JJ  =  1 

200  READ  (5  ,106)MDEG,(  C(l),  1*1,  6) 

I F ( MDEG.GT • 6 )  GO  TO  48 

WR  I  TE  (  fet  145 )  <  l,C( I ), I  =  l,MDEG) 

IF  ( MDEG  -  6)  49,  49,  48 

48  READ  (5  ,104)(CII),  I  *  7,  MDEG) 

WR I TE ( 6» 145 )  (  I,  CII ), 1=1, MDEG) 

106  FORMA  T (  I 10.6EIO.O) 

49  IF  (  I GA ( J J  ) -  1  )  50,  51,  50 

51  DO  56  1=  1,  MDEG 

56  ANSI  I  )  *  CII) 

[ADEG  =  MDEG 

IF  (ISDEG  -  MDEG)  52,  52,  53 

53  INDEG  *  ISOEG 
GO  TO  68 

52  INDEG  *  MDEG 
GO  TO  68 

50  READ  (5  ,106 )NDEG,  (Dili,  I  *1,  6) 

IF ( NDEG.GT • 6 )  GO  TO  54 
WRITEI6.145)  I  I, D(I ), |*1,N0EG) 

IF  (NOEG  -  6)  55,  55,  54 

54  READ  (5  , 104 ) ( D( I ) ,  1*7,  NDEG ) 

WRITEI6.145)  (  I, DU),  1*1, NOEG) 

55  I ADEG  *  NDEG  ♦  MDEG  -l 

CALL  POLMPY  ( C, MDEG, D, NDEG, ANS ) 

IGA  (JJ)  *  IGA(JJ)  -1 
IF  (IGA(JJ)  -l)  65,  65,  64 

64  DO  60  I  *  l,  I ADEG 
60  C( I )  *  ANSI  I  ) 

MDEG  *  I ADEG 
GO  TO  50 

65  IF  (ISOEG  -  I ADEG )  66,66,67 

66  INDEG  *  I ADEG 
GO  TO  68 

67  INDEG  *  ISOEG 

68  CALL  POLADD  ( SAV E , I SDEG , ANS , I  A OE G , ER A ) 

145  FORMAT (  10X,2HC(,  I  3, 2H ) * , F  20. 10 ) 

WR  l  TE  (  6  ,  l) 

1  FORMAT (  IH 1 / / ) 

6800  IF  (ERA( INDEG) )6803,  6802,  6803 
6802  INDEG  =  INOEG  -  l 

IF  (INOEG)  6801,  6801,  680C 


67 

68 
145 

1 

6800 

6802 

6801 

6803 

LOO 

r"- 

oo 

LO 


INDEG  *  1 
JJ  =  JJ  ♦  l 
00  70  I  =  1,  INDEG 
SAVE! I  )  =  ERA (  I  ) 

ISDEG  =  INDEG 
IA  *  IA  -l 

IF  ( IA)  201,  201 ,  200 
IF  ( FORK  l )  202,  202,  203 
SAVE  NUMERATOR. 

00  220  I  =  1,  ISDEG 


i 


! 
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220  AU)  =  SAVE  (  I  )  Ini 

IDA  *  ISOEG  M2 

IA  *  IB  U3 

FORK l  »  1.0  lri4 

00  230  1  =  1,  IA  IMS 

230  IGA ( I )  s  IGB ( I )  loo 

C  START  DENOMINATOR  1  1 1 

GO  TO  300  lb* 

C  SAVE  DENOM.  U9 

203  DO  240  I  *  1,  ISOEG  l 40 

240  HIT)  =  SAVE!  1  )  Ml 

IDB  *  ISOEG  1 42 

WRITE  (6  ,  IC9)  I'M 

109  FORMAT  (  10X.41HC0EFF I  Cl  ENTS  ARE  GIVEN  IN  ASCENDING  ORDER///// )  M4 

339  IF  (A(IDA))  340,  341,  340  MS 

341  IDA  *  IOA  -  1  1  it> 

IF  I  IDA)  345,  345,  339  1  i 1 

345  WRITE  (6  ,120)  Mb 

120  FORMAT  (  1H0, 10X, 2OHPOLYN0MI AL  A  IS  ZERO//)  199 

F0RK2  *  1.0  200 

GO  TO  410  201 

340  IF  (IDA  -  2)  346,  347,  335  202 

346  WRITEI6, 121)A(1)  203 

STAR*ST  AR 1  2u4 

PRNT*.TRUE.  205 

ANUMB1»A(1)  206 

ANUMB2«0.0  207 

CALL  PL0TER(SAVE1,SAVE2, A NUMB 1 , A NUMB 2 , POINT .XA.XP.MI.PRNT.iTAR, I  II  20* 

*,N0)  209 

121  FORMAT  (1H0,10X,?8HP0LYN0MIAL  A  IS  A  CONSTANT  * , 1 P IE  16. 7// )  210 

GO  TO  410  211 

347  ROOT  «  -  All)  /  AI2)  212 

WRITE  (6  « 133 ) A( 1),  AI2)  213 

133  FORMAT  (10X,21HTHE  COEFFICIENTS  OF  A/1P2E20.7)  214 

WRITE  (6  , 122 )ROOT  215 

ST  AR*ST  AR 1  216 

PRNT=.TRUE.  217 

ANUMB1-R00T  21H 

ANUMB2=0.0  219 

CALL  PLOTERI SAVE  1, SAVE2, ANUMB1 , ANUMB2 , POINT ,XA,XB,MI,PRNT,STAR,I 1 1  220 

♦.NO)  221 

122  FORMAT  I IHO, 10X, 23HR00T  OF  POLYNOMIAL  A  I S , l PI E16. 7// )  222 

GO  TO  410  223 

C  WRITE  POLYS  224 

335  ID1A  = I  DA  -l  225 

WRITE  (6,107)ID1A,(A(I),I«1,I0A)  226 

K  »  IDA  227 

DO  BOO  I  *  1, IOA  228 

ASCII  *  AIK)  229 

800  K  *  K-l  230 

IDP2A* IOA  *2  231 

ID2A*  2  ♦ I D IA  2 32 

CALL  MULLER  (AS,  ID1 A.ROOTR.ROOTI )  233 

DO  805  I  =  l.IDIA  234 

SAM  *  100.  ♦  AMAXKABSIROOTRI  I  ))  .ABSIROOTI  (I  )  )  )  235 

IF  (SAM  *  ABSIROOTRI I ) ).EQ.  SAM)  RCOTR 1 1 ) «  0.0  236 

IF  (SAM  ♦  ABSIROOTI 1 1 ) ) •  EQ.  SAM)  RCCT I  III*  0.0  237 

805  CONTINUE  23B 

400  WRITE  (6,111)  (RCOTR I  I) .ROOT  I  I  I)  ,  I  =  1 , 1 Dl A )  2 39 

CALL  ERCHEK(ROOTI.IOIA)  240 

64 
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< 


O  PRNT-. TRUE.  2*1 

STAR-STAR  1  2*2 

00  2  I I l-l , ID1A  2*3 

ANUMB1-R00TR( III  )  2*4 

ANUMB2-ABSIR00TI (111))  2*5 

CALL  PLOTER! SAVE l, SAVE2 ,  ANUMB1 ,ANUMB£  » POINT ( XA ,XB » Ml »PRNT , STAR, I II  2*6 

*,N0)  2*7 

2  CONTINUE  2*8 

*10  IF(BUDB))  *11,  *12,  *11  2*9 

*12  IDB  *  IOB  -  1  250 

IF  (IOB)  **5,  **5,  *10  251 

**5  WRITE  (6  ,123)  ?52 

123  FORHAT  ( IHO, 10X, 20HPOLYN0MI AL  B  IS  ZERO//)  253 

IF  (F0RK2) 12, *50, 12  25* 

*50  F0RK3  *  1.0  255 

GO  TO  698  256 

*11  IF  (IDB  -  2)  *51,  *52,  *99  25/ 

*51  WRITE  (6  , 12* ) B( IOB )  258 

STAR«STAR2  259 

PRNT*.) RUE.  260 

ANUMB1-B( IDB)  261 

ANUMB2-0.0  262 

CALL  PLOTER ( SAVEl ,SAVE2 , ANUHB1 ,ANUHB2 , POINT ,XA,XRVMI , PRNT, STAR , III  263 

*,N0)  26* 

PRNT-.FALSE.  265 

STAR-STAR3  266 

12*  FORMAT  ( IHO, 10X, 28HP0LYN0MI AL  B  IS  A  CONSTANT  », 1PIE 16. 7// )  267 

GO  TO  698  268 

*52  ROOT  -  -BID  /  B(2)  269 

WRITE  16  «13*)B( 1),  B  ( 2  )  270 

13*  FORMAT  ( 10X» 21HTHE  COEFFICIENTS  OF  B/IP2E20.7)  271 

WRITE  16  , 125 )ROOT  272 

STAR-STAR2  273 

PRNT-.TRUE.  27* 

ANUMBl-ROOT  275 

ANUMB2-0.0  276 

CALL  PLOTER (SAVEI,SAVE2,ANUMB1,ANUMB2, POINT, XA,XB,MIV PRNT, STAR, I II  2  77 

*,N0)  278 

PRNT*. FALSE.  279 

STAR-STAR3  200 

125  FORMAT  (  IHO, 10X, 23HR00T  OF  POLYNOMIAL  B  I S, 1 P1E16.7// )  281 

GO  TO  698  282 


I 


107  FORMAT  (10X,*2HTHE  COEFFICIENTS  OF  POLYNOMIAL  A  (ORDER  »  13, 1H)/  (  283 

11P6E20. 7) )  28* 

*99  1010  •  IOB  -l  285 

WRITE  (6, 108)ID1B,(B(I), 1-1,108)  286 

108  FORMAT  (////10X,*2HTHE  COEFFICIENTS  OF  POLYNOMIAL  B  (ORDER  *  13, 1H  287 

1)/  (1P6E20.7) )  288 

K  «  IOB  289 

00  801  I  -  1,108  290 

BS(I)  «  B(K)  291 

801  K  -  K-l  292 

I0P2B-  IOB  *  2  293 

1028  -2*  1018  29* 

CALL  MULLER  ( BS, ID1B,R00TR,R00TI )  295 

00  806  I  -  1,1018  296 

SAM  «  100.  *  AMAX1(ABS( ROOTRI I )) ,A8S(R00TI ( I ) ) )  297 

IF  (SAM  ♦  ABS( ROOTRI I ) ) «E0.  SAM)  ROOTR(I)*  0.0  298 

IF  (SAM  ♦  ABSIROOTl ( I ) ) • EQ.  SAM)  ROOTI ( I )-  0.0  299 

806  CONTINUE  ,,  300 


58767 


500  WRITE  ( 6, 1 12 ) ( R00TRI I ), ROOT  HI ) , I  *  1.ID1B)  301 

CALL  ERCHEMROOTI,  ID1B)  3u2 

STAR-STAR2  30 3 

00  3  111*1, ID1B  304 

ANUMB1-R00TR1 III )  305 

ANUMB2*ABS(R00TI ( III) )  306 

CALL  PL0TER(SAVE1,SAVE2, ANUMB1 ,ANUHB2,P0INT,XA,XR, MI «PRNT«ST AR» III  307 

♦  »N0)  30(3 

3  CONTINUE  309 

PANT*. FALSE,  310 

STAR«STAR3  3 1 1 

111  FORMAT  (lHOtllX* 14HTHE  ROOTS  OF  A/  (1PIE20.7.6H  41  , IPIE14. 7, 1P1  312 

1E20.7.6H  ♦!  .1P1E14.7, 1PIE20.7.6H  +  1  ,1P1E14.7))  313 

112  FORMAT  ( lHO* UX, 14HTHE  ROOTS  OF  B/  (1P1620.7,6H  ♦!  , IP1E14. 7, 1P1  314 

1E20.7.6H  41  , 1P1E14.7, 1PIE20.7,6H  41  .1P1E14.7I)  315 

698  IF  I  FORK 2 ) 12.699. 12  316 

699  IF  (F0RK3) 12.6991. 12  317 

6991  WRITE  (6  ,102)IPR0B  318 

MSHEET  *  5  319 

C  START  K  CALCULATIONS  320 

IF  (N)  702.702,533  321 

533  00  550  I*  1.  N  322 

00  541  J*  1.  IDA  323 

541  ATM  J)  *  XIII  *  A(  J )  324 

C  COMPUTE  ROOTS  Of  K  *  A  ♦  B  325 

IOC*  MAXOI  IDA.  I DR I  326 

CALL  POLADD  ( ATK, IDA.B, I  OB,  CM  327 

IOS  «  IOC  328 

554  IF  (CM  IOS)  1555,  557,  555  329 

557  IDS  *  IDS  *  1  330 

IF  (IDS)  558,558,  554  331 

558  WRITE  (6  ,129)X(I)  332 

129  FORMAT  1 1H0, 10X,35HP0LYN0M!AL  K*A  4  B  IS  ZERO  FOR  K  * , 1PI£ 1 6. 7// )  333 

GO  TO  550  334 

555  IF  (IDS  -  2)  559,  560,  561  335 

559  WRITE  (6  , 130JCKC IDS) •  X<!)  336 

130  FORMAT  1 1H0, 10X, 35HP0LYN0MI AL  K*A  4  B  IS  A  CONSTANT  *  , IP1E15.7.10  337 

IH  FOR  K  *  , 1P1E14. 7//)  338 

GO  TO  550  339 

560  ROOT  «  -CX(l)  /  CK( 2)  340 

WRITE  (6  ,  131 IROOT,  XU)  341 

131  FORMAT  ( 1H0, 10X, 18HR00T  OF  K*A  ♦  B  ■  , 1P1E15.7, 10H  FOR  K  *  , IP1EI  342 

14.7//)  343 

GO  TO  550  344 

561  K  -  IDS  345 

00  803  J  -  1, IOS  346 

CKS( J)  *  CMK )  347 

803  K  -  K  -  1  348 

I01C  «  IOS  -  l  349 

I0P2C  *  IDS  ♦  2  350 

I02C  *24  IDIC  351 

CALL  MULLER  (CXS.IOIC.ROOTR.ROOTI )  352 

00  807  «)  «  1.I01C  353 

SAM  -  100.  6  AMAXK  A8S(R00TR(  J))  ,  ABSC  ROOT  I  ( J ) ) )  354 

IF  (SAM  4  A6SCROOTRI J) I.EQ.SAM)  ROOTR(J)  -  0.0  355 

IF  (SAM  4  ABSIROOTK J ) I.EQ.SAM)  ROOTI(J)  ■  0.0  356 

807  CONTINUE  357 

WRITE  (6,808) IDIC,X(1)«(CK(J),J*1,I0S)  358 

808  FORMAT  ( ///10X.48MTHE  COEFFICIENTS  OF  POLYNOMIAL  K*A  4  B  (OROER  -  359 

113, 7H)  K  •  1P1E16.7/UP6E20.7))  360 


O  545  WRITE  ( 6» 115) (ROOTR(J) » ROOT  I (J)«J*1,1DLC)  <61 

CALL  ERCHbKIROOT I « I 0 1 C )  36? 

CALL  SAVER! ROOTR , ROOT I» 10 1C* SAVE  l,SAVE2,JZ0,Kl )  <6  3 

115  FORMAT  ( IHO,9X,16HROOTS  OF  K*A  4  B/ ( IP1F20. 7 , 6H  ♦  I  , 1PIF14. 7, IP l  364 

1E20.7,6H  4  I  .IP1E14.7, 1P1E20.7,6H  ♦  I  ,1P1E14.7)»  365 

5452  MSHEET  *  MSHEET  -  1  366 

IF  (MSHEET)  546,  546,  550  367 

546  WRITE  (6  ,102)!PROB  36R 

MSHEET  *  5  36) 

550  CONTINUE  37u 

GO  TO  12  371 

702  DO  705  J  *  1,  IDA  3/2 

705  ATM  J  )  =  Y  *  A I J  )  373 

C  COMPUTE  ROOTS  OF  K  *  A  4  B  374 

IOC*  MAXOUDA,  IDB)  375 

CALL  POLAOO  ( ATK , IDA,  B, I DB,CK )  376 

IDS  *  IOC  377 

754  IF  (CM  IDS  I  )  755,  757,  755  378 

757  IDS  *  IDS  -  1  37) 

IF  (IDS)  758,  758,  754  <80 

750  WRITE  (6  , 129)Y  381 

GO  TO  711  382 

755  IF  (IDS  -  2)  759,  760,  761  383 

759  WRITE  (6  , 130)CK( IDS) ,  Y  384 

GO  TO  711  385 

760  ROOT  *  -CM  1 )  /  CM2)  386 

WRITE  (6  ,131 )ROOT ,  Y  387 

GO  TO  711  388 

761  K  *  IDS  389 

DO  804  I  *  1,  IDS  390 

CK SCI)  *  CK ( K  )  391 

804  K  *  K  -l  39? 

I01C  *  IDS  -  l  393 

IDP2C  *  IDS  *  2  394 

ID2C  *  2  *  I01C  395 

CALL  MULLER  ( CKS, ID1C , ROOTR , ROOT  I  )  396 

DO  009  I  =  l.IDIC  397 

SAM  *  100,  *  AMAX1 ( ABS ( ROOTR ( I )) ,ABS ( ROOT  1(1)))  398 

IF  (SAM  4  ABS ( ROOTR ( I ) ) • E Q •  SAM)  ROOTR(I)*  0.0  399 

IF  (SAM  4  ABS (ROOT I ( I ) ) .EQ.  SAM)  ROOT I ( I )*  0.0  400 

809  CONTINUE  401 

WRITE  ( 6, 808) I01C,Y,(CK( I ), 1*1, IDS)  402 

WRITE  (6,115)  (ROOTR(J) , ROOT  I ( J) ,J*1,ID1C)  4  3 

CALL  ERCHEM  ROOT  I  * ID1C )  404 

CALL  SAVER ( ROOTR , ROOT I ,ID1C,SAVEI,SAVE2,JZ0,K1)  405 

711  Y  »  V  4  DY  406 

IF  (Y  -  YT)  712,712,12  407 

712  MSHEET  «  MSHEET  -  1  400 

IF  (MSHEET)  713,  713,  70?  409 

713  WRITE  (6  , 102 ) IPROB  410 

MSHEET  *  5  411 

GO  TO  702  412 

0&016  CALL  PL0TER(SAVEI,SAVE2,ANUMB1 ,ANUMB2, POINT, XA,XP, Ml, PRNT, STAR, i II  413 

CD  *,N0)  414 

Is-  READ(  5, 22  )L  415 

°0  22  FORMAT! I 10 )  416 

10  IF(L.EQ.O)  GO  TO  1017  417 

CALL  EXPAND! L , XA , XB )  418 

1017  RETURN  419 

END  420 

67 
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c 

c 

c 

c 

SUBROUTINE  POLMPY  (AtNtB,M,C) 
DIMENSION  MllfBIll.CIl) 

K  «  M*N 
00  5  I»1,K 

5  cm  *  0.0 

00  10  I*1,N 
L  -  1-1 
00  10  J*1,M 
L  *  L>1 

10  C(L )  -  C(L>*A(I)*B(J) 

RETURN 

END 

C 

C 

C 

C 

SUBROUTINE  P0LA00  (A,N,B,M,C) 
DIMENSION  Alll.BIlltCm 
IF  IN-M)  1,1,2 

1  NK  »  N 
00  TO  5 

2  NK  *  M 

5  00  10  I«1,NK 

10  cm  «  Alikin 

NK  *  NK*I 
IF  IN-MI  11,23,15 

11  DO  20  I*NK«M 
20  CIII  •BID 
25  RETURN 

15  00  30  I«NK«N 

30  cm  -  All) 

RETURN 

END 

C 

C 

C 

C 

SUBROUTINE  MULLERICOE ,N1,R00TR,R00TI I 
DIMENSION  COE 1 1 1 ,R00TR 1 1 ) , ROOTI 1 1 1 
N2*N1+1 
N4»0 
I-NU1 

19  IF(COEII) 19,7,9 
7  N4*N4+1 
R00TRIN4 1*0* 

KCOTKN4I-0. 

1*1-1 

IFIN4-N1I 19,37,19 
9  CONTINUE 
10  AXR-0.8 
AXI-O. 


L*1 

N3»l 

ALP1R-AXR 

ALP1I-AXI 


421 

422 

423 

424 

425 

426 

427 

428 

429 

430 

431 

432 

433 

434 

435 

4  36 

437 

4  3H 

439 

440 

441 

442 

443 

444 

445 

446 

447 

448 

449 

450 

451 

452 

453 

454 

455 

456 

457 

458 

459 

460 

♦6 1 

HPRS  462 

463 
HPRS  464 
HPRS  465 
HPRS  466 
HPRS  467 
HPRS  468 
HPRS  469 
HPRS  470 
HPRS  471 
HPRS  472 
HPRS  473 
HPRS  474 
HPRS  475 
HPRS  476 
HPRS  477 
HPRS  478 
HPRS  479 
HPRS  480 
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G0T099 

11  BET1R-TEMR 
BET1I-TEMI 
AXR*0.85 
ALP2R-AXR 
ALP2I-AXI 
M«2 

G0T099 

12  BET2R-TEMR 
BET2I-TEHI 
AXR-0.9 
ALP3R-AXR 
ALP3I-AX1 
M«3 

G0T099 

13  BET3R-TEMR 
BET3I-TEMI 

14  TEl*ALPlR-Al.P3R 
TE2* ALP l I -ALP 3 I 
TE5-ALP3R-ALP2R 
TE6-ALP3I-ALP2I 
TEH-TE54TE54TE64TE6 
TE3«| TEl*TE5*TE2*TE6) /TEM 
TE4*I TE2*TE5-TEl*TE6) /TEM 
TE7-TE3+1. 

TE9-TE34TE3-TE44TE4 
TE10«2.*TE3*TE4 
DE15«TE7*BET3R-TE4*BET3 I 
0E16-TE7*BET3DTE44BCT3R 

TE 1 1-TE3*BET2R~?E4*BE  T2 1 ♦BET1R-OE15 
TE12-TE3*BET2  DTE4*BE  Y2R+BET1I -OE 16 
TE?«TE9-1. 

TE1»TE9*BET2R-TE 10*BET2 I 

TE2-TE9*BET2I+TE10*BET2R 

TE13«TE1-BETIR-TE7*BET3R*TEL0*BET31 

TE14«TE2-BET1I-TE7*BET3I-TE10*BET3R 

TE15«DE15*TE3-0E16*TE4 

TE16»0E15*TE44DE 16*TE3 

TE1»TE13*TEL3-TEL4*TE14-4.*(TE11*TE15-TE12*TE16I 
TE2*2.4TE13*TE 14-4.61  TE 12*TE15+TE1 l*TE16) 

TEM  *  SORT  I  TEI*TE1«-TE2*TE2) 

IF( TE1)  113,113,112 
113  TEA  ■  SORT  1.5  *  ITEM  -  TE1D 
TE3».5*TE2/TE4 
GO  TO  111 

112  TE3  =  SORT  1.5  6  ITEM  ♦  TED) 

IFITE2D10, 200,200 
DO  TE3--TE3 
200  TE4*.5*TE2/TE3 
111  TE7«TE13tTE3 
TEB»TE144TEA 
TE9*TEI 3-TE  3 
TE10«TE14-TE4 
TE1-2.4TE15 
TE2«2.*TF16 

IFITE74TE74TE8*TE8-TE96TE9-TE104TE1C) 204,204,205 

204  TE7-TE9 
TER*TE10 

205  TEM«TE7*TF7*TE8*TE8 
TE3»CTE1*TE74TE24TE8)/TEM 


HPRS 

481 

HPRS 

482 

HPRS 

483 

HPRS 

484 

HPRS 

485 

HPRS 

486 

HPRS 

487 

HPRS 

488 

HPRS 

489 

HPRS 

<90 

HPKS 

4)1 

HPRS 

<  72 

HPRS 

*93 

HPKS 

494 

HPRS 

495 

HPRS 

496 

HPRS 

497 

HPRS 

498 

HPRS 

499 

HPRS 

SCO 

HPRS 

501 

HPRS 

502 

HPRS 

503 

Hf  KS 

504 

HP.7S 

505 

HPRS 

506 

HPRS 

507 

HPRS 

508 

HPRS 

509 

HPRS 

510 

HPRS 

511 

HPRS 

512 

HPRS 

513 

HPKS 

514 

HPRS 

515 

HPRS 

516 

HPRS 

517 

HPRS 

518 

HPRS 

519 

HPRS 

520 

HPRS 

521 

HPRS 

522 

HPRS 

523 

HPRS 

524 

HPRS 

525 

HPRS 

526 

HPRS 

527 

HPRS 

528 

HPRS 

529 

HPRS 

530 

HPRS 

531 

HPRS 

532 

HPRS 

533 

HPRS 

534 

HPRS 

535 

HPRS 

536 

537 

538 

HPRS 

539 

HPRS 

540 
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t 


( 


( 


( 


K 


TE4«( TE24TE7-TE14TE8) /TEH  HPRS  541 

AXR«ALP3R4TE3*TFS~<E4*TE6  HPRS  542 

AXI*ALP3I+TE3*TE6*TE4*TE5  HPRS  043 

ALP4R-AXR  HPRS  544 

ALP4I*AXI  HPRS  545 

M«4  HPRS  54ft 

60  TO  99  HPRS  547 

15  N6*l  54 h 

38  IF  IABS  (HELL)  ♦  ABS  (BELL)  -  l.E-20)  18,18,16  549 

16  TE7  •  ABS  ( ALP3R  -  AXR)  ♦  ABS  ( ALP3I  -  AXI)  550 

IF  (TE7  /  (ABS  (AXR )  ♦  ABS  UXD)  -  l.E-7)  1 8, 18, 17  551 

17  N3*N3+i  552 

ALPIR-ALP2R  553 

ALPI1-ALP2I  554 

ALP2R-ALP3R  555 

ALP2UALP31  556 

ALP3R*AL('4R  55  7 

ALP3I-ALP4I  658 

8ET1R-6ET2R  557 

BETII-6ET2I  560 

BET2R-BET3R  561 

BET2I-BET3I  562 

BET3R*TEHR  563 

BET31»TEMI  564 

IF(M3-100) 14, 18, 18  565 

18  N4-N4+1  56ft 

ROOTR ( N4 ) * ALP4R  567 

R00TKN41-ALP4I  568 

H3*0  569 

41  IF(N4-Ni ) 30,37,37  570 

37  RETURN  571 

30  IF  (ABS  (R00TKN4) )  -  l.E-5)  10, 10,31  572 

31  GO  T0( 32, 10), L  573 

32  AXR-ALP1R  574 

AXI— ALP1I  575 

ALP1I—  ALPil  576 

H»5  577 

60  TO  99  578 

33  BET1R-TENR  579 

BETU-TEMI  580 

AXR-ALP2R  581 

AXI— ALP2I  582 

ALR2I— ALP2I  583 

H*4  584 

GO  TO  99  585 

34  BET2R-TEHR  586 

BET2UTEM1  587 

AXR-ALP3R  588 

AXI— ALP3I  589 

ALP3I— ALP3I  590 

L«2  591 

M*3  592 

99  TEHR«COE( 1 I  593 

TENI-0.0  594 

001001*1, N1  595 

TEt*TENR*A;<R-TEHI*AXI  596 

TEH1»TEHI*aKR4TEHR*AXI  597 

100  TERR*  TE14C0EII41)  598 

MELL*TEHR  599 

8ELL-TEHI  600 


70 
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c 

c 

c 

c 


c 

c 

c 

c 


42  IF! N4)l 02 *103*102 

102  001011*1* N4 
TEM1*AXR-R00TR(I ) 
TEM2*AXI-R00TI(  I ) 
TEI*TEM14TEM14TEM2*TEM2 
TE2*(TEMR*TEM1«TEMI*TEM2 ) /TE1 
TEMI*(TEMI*TEHl-TEMR4TEM2)/TEl 

101  TEMR-TE2 

103  GO  T0!ll*12*13*l5,33,34),M 
END 


SUBROUTINE  ERCHEK(X*1) 

DIMENSION  XtlOO) 

DATA  ERllLIM/O* IE-8/ 

00  10  J»1*I 

10  IFIABSIXI J)  ULT.ERRLIM)  X(J)*0*0 
RETURN 
END 


SUBROUTINE  EXPAND! I LOVE »  X,Y I 

DIMENSION  X2( 2000)*Y2(2000)  *X!2000) *Y(?000) 
DIMENSION  Xl(  2000), YU 20001  *2!  1251*21(87) 

INTEGER  SPOT! 130*50) *DASH,Q, BLANK, POINT! 130, 100) 
OATA  Q/1HI/,  OASM/lH-/,BLANK/lM  / 

OATA  APLUS/O*  O/.BPLUS/O*  0/*AMI NUS/0*0/ ,BMI NUS/0*0/ 

COMMON  / INF04/  0, BLANK* DASH 

COMMON  /BOB/POINT 

COMMON  /PAYNE/AM, AP,BM,BP 


COMMON  /PJ/X1 » Yl 
EQUIVALENCE (XI! 1 )*Z!1I),!Y1(1),21(1)) 
EQUIVALENCE  ( SPOT! 1* 1 )*P0INT( 1 *1 ) ) * ! I l* 12 
EQUIVALENCE! AP*APLUS)*( AM* AMI NUS I  *  !  BM*BMIl 
)0  40  J*l* ILOVE 


*  i 

BMINUS)* (BP.BPLUS) 


30  50  M*  1,2000 
!1(M)*0#0 
ri(M)-o.o 
30  51  M*l,  130 
30  51  N*l,  50 
SPOT!M,N)*BLANK 
30  52  M*l,  130 
SPOTIM*  1  )*OASH 
SP0T!M,50)*DASH 
30  53  M*l,50 
SPOT!  1,M)*Q 
SPOT!  126,M)*0 
DO  54  M*  1,2000 
K2(M)*0«0 
72(M)*0*0 

READ! 5,10)  OMEGA  *  ENCRMT , S IGMA  *  DE  LTA 
DMEGA-ABS IOMEGA) 

FORMAT (4F 10*0) 

A-ENCRMT 


B-OMEGA 


71 


601 

60? 

603 

604 

605 

606 

60  7 
608 
60S 
610 
611 
612 

61  J 
614 
6 1  5 
616 

617 

618 

619 

620 
621 
622 

623 

624 
62b 
626 

627 

628 

629 

630 

631 

632 

633 

634 

635 

636 

637 

638 

639 

640 

641 

642 

643 

644 

645 

646 

647 

648 

649 

650 

651 

652 

653 

654 

655 

656 

657 

658 

659 

660 


58773 


/ 
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c 


( 


( 


r 


B-B'PERCNT  661 

APLUS-B40NE6A  662 

AHINUS-ONECA-B  663 

C-APLUS  664 

O-ANINUS  66b 

11*1  666 

00  41  L-l»2000  667 

IF(Y<L)«LT«0)G0  TO  41  668 

1FCYIL )«GT«C)GO  TO  41  669 

23  XUUI-XILl  670 

YlUl)-Y(L)  671 

11- Il+l  672 

41  CONTINUE  673 

A-OELTA  676 

901  PERCNT-A*O.Ol  67b 

B- SIGMA  676 

903  6-B*PERCNT  677 

BPLUS-SIGMA+B  678 

BNINUS-SIGNA-B  679 

C-BPLUS  680 

O-BMINUS  681 

12*1  68? 

910  00  42  L-l* 2000  683 

IF(C«GT«0.0)  GO  TO  60  684 

IF(XHL).GT.D)  GO  TO  42  685 

IFIX1CU.LT. C)  GO  TO  42  686 

GO  TO  25  687 

*0  IFCXl(L).GT.C)  GO  TO  42  688 

IF(XUL).LT.D)  GO  TO  42  689 

25  X2(I2)-X1(L)  690 

V2(I2)-Y1(L)  691 

12- 1241  692 

42  CONTINUE  693 

CALL  SPLIT(X2tY2t SPOT«APLUStAM 1NUS» 6PLUS tBMINUS)  694 

40  CONTINUE  695 

RETURN  696 

END  697 

C  698 

C  699 

C  700 

C  701 

SUBROUTINE  SPLITfX»Yv SPOTtAPLUS* AMINUS«BPLUSVBMNUS)  702 

OIMENSION  Xl(2000)« Y1I2000)  703 

OINENSION  X(  2000 ) »Yt 2000) «2( 125) • Z1C87)  704 

INTEGER  SP0T(130t50)  705 

CONNON  /PJ/Xl* VI  706 

EQUIVALENCECX1I1 ) «2(  1))»(Y1(1)«21(1))  707 

A-APLUS  708 

B-AMINUS  709 

C-A-B  710 

O-BPLUS  711 

E-BNINUS  712 

G-ABSIOI  713 

H-ABS(E)  714 

F-G-H  715 

PELTA-C/124.0  716 

0IFF-F/50.0  717 

00  11  J-U124  718 

2(J)-A  719 

II  A-A-OELTA  720 


72 
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o 


Z( 12b  1*8  721 

IF!O.LT.O.O)OIFF«- OIFF  72? 

00  12  J*l*49  723 

Z1(J)*D  724 

12  OO-DIFF  7 2b 

ZI(50)*E  726 

CAIL  BRAKUP  (  Z  •  X«  Y»  SPOTt  Z  1 )  72  7 

RETURN  728 

ENO  729 

C  730 

C  711 

C  73.? 

C  713 

SUBROUTINE  BRAKUP< YY*Xf V, SPOT* XX )  734 

LOGICAL  SKIP*  SKI  PI  73b 

INTEGER  SPOT! 130»b0) »  XPT,YPT,STAR  716 

DATA  STAR/lH*/  737 

OIMENSI ON  X(2000>,Y(2000>,YY! l25)fXX(87)  71M 

COMMON  /PAYNE/AHf AP»BM*BP  719 

EQUIVALENCE! APfAPLUS) » ( AM»AMINUS ) « ( BMvBMINUS ) « I  DP VB PLUS )  740 

SKIP««FALSE.  741 

SKIP1-. FALSE*  74? 

00  43  J«1 * 2000  743 

00  43  I*l« 125  744 

IFIVI J).NE*0*0)  GO  TO  30  74b 

IFC X( Jl.NE.O.O)  GO  TO  30  746 

L*J*1  74  7 

M*J*6  748 

00  60  N«L.M  749 

IFIN.GT.2000)  GO  TO  60  7S0 

IF(Y(N)*NE*0*0)  GO  TO  30  7bl 

IF(X(N)«NE*0*0)  GO  TO  30  7b2 

60  CONTINUE  7b3 

GO  TO  40  754 

30  IF(SKIP)  GO  TO  20  Tbb 

1FCYC J).LT.YY(  III  GO  TO  20  7b6 

YPT-I  7b 7 

SKIP«*TRUE*  7brt 

20  IFC  SKIP1 )G0  TO  48  7b9 

IFIWGT.50)  GO  TO  43  760 

IF(8PLUS*GT*0*0)  GO  TO  10  761 

IFtXU)*GT*XXim  GO  TO  43  762 

GO  TO  11  763 

10  mxi j).lt.xx<  in  go  to  43  764 

11  I F( BPLUS*LT*0*0)  XPT-I  76b 

IFIBPLUS.GE.O.)  XPT*5l- 1  766 

SKIP1-.TRUE.  767 

IF| *NOT*SK IP )G0  TO  43  76U 

50  1*125  769 

SPOTI VPT* XPTI*STAR  7/0 

SKIP-*FALSE*  771 

SKIPI*«FALSE*  772 

48  IFISKIPIGO  TO  50  773 

43  CONTINUE  774 

40  CALL  RITEITISPOT)  77b 

RETURN  776 

ENO  777 

778 

779 

780 
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C  781 

SUBROUTINE  RITEITISPOT)  7 82 

COMMON  /PAYNE/AM, AP,BM*BP  783 

EQUIVALENCE! APvAPLUS)t<A«(AMf  *>)  ,(BH,BMINUS) , (BP.BPLUS)  784 

INTEGER  SPOT! 130*50)  785 

IF(BP.LT.BM)  GO  TO  40  78o 

A«BP  787 

B*BM  788 

MRI TE16, 1)  789 

1  FORMAT! 1HI,40X,43HC0MPLEX  FREQUENCY  PLANEtRIGHT  HANO  QUADRANT  )  790 

GO  TO  30  791 

40  A*BM  79? 

B»BP  793 

MRITE!6,2)  794 

2  FORMAT! 1H1,40X,43HC0MPLEX  FREQUENCY  PLANE,  LEFT  HAND  QUADRANT  )  795 

30  WRITE(6,12)AP,AM  796 

12  FORMAT! 1|4X«15H< - J-0MEGA,/1X,F8.2, 112X, F8.2)  797 

WRITE! 69 14)B  798 

14  FORMAT ( 124X,5HSIGMA,/126X,1HI « /126X,1HV,/121X,F8«2)  799 

WRITE16, 11)  SPOT  800 

11  FORMAT! IX, 130AI)  B01 

WRITE(6« 15) A  802 

15  FORMAT!60X*27HLINEAR  EXPAND  PLOT! RAD/SEC >  »34XVF8«^)  803 

RETURN  804 

END  805 

C  806 

C  807 

SUBROUTINE  SAVER IROOTR* ROOT I • 1 01C«SAVE1, SAVE2,  JZ0»K1 )  808 

DIMENSION  SAVE1! 100, 100 1 «SAVE2 (100*100) tROOTR! 1001. ROOT  I (100)  809 

IF! Kl )30,9, 10  810 

9  Kl»l  811 

10  I01M01C  812 

ID1C-ID1C+!K1>1)  813 

IFUOIC.GE.IOOIGO  TO  30  814 

50  00  40  IZAP»K1,ID1C  815 

IZAMZAP-(K1  -1)  816 

SAVE1UZ0,IZAP)-R00TR!IZA)  817 

SAVE2(JZ0, IZAP)*ROOTI!IZA)  818 

40  CONTINUE  819 

K1**ID1C+1  820 

GO  TO  20  821 

30  JZ0*JZ0+1  8 22 

IDIC*ID1  823 

GO  TO  9  824 

20  RETURN  825 

END  826 

C  827 

C  828 

C  829 

C  830 

SUBROUTINE  PL0T2R  ( SAVE1  •  SAVE2  *  ANUMB1 ,  ANUMB7.  *P0 1  NT  ,  XA  *  XB  *  MI  *  PRNT  *  831 

*ST*R,NBB,NO)  832 

OIMENSION  SAVE1( 100,100) ,SAVE2(100, 100), XA(200G) ,XB(2000)  833 

INTEGER  POINT! 130,100) tOASH, BLANK, STAR  834 

LOGICAL  PRNT  835 

IF! PRNT)  GO  TO  43  836 

00  40  NAB»1,100  837 

00  40  NBB>1,100  838 

IF! SAVE1(NAB«NBB )) 41 ,42, 41  839 

42  IF(SAVE2(NAB,NBB) 141,40,41  840 
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1 


o 


41 


43 


50 

45 

40 

70 


C 

C 

C 

C 


500 


SO 


c 

c 

c 

c 


21 


40 


31 


ANUMB1-SAVE1(NAB*NBB) 

ANUMB2-S AVE2 * NAB*NBB 1 
ANUMB2-ABS*  ANUM62 1 
AZZ-ANUMB2 
AZZ-AZZ*100000.0 
NZZ-AZZ 

IFINZZ.EQ.O)  GO  TO  45 
NBB-NBB+1 

CALL  EXCUTEt ANUMB 1* ANUMB 2, POINT* XA*XB* MI , NO, STAR) 

IF(PRNT) RETURN 

CONTINUE 

CALL  MR1TIT*  XA«X6) 

CALL  PREPARE  POINT  I 

RETURN 

END 


SUBROUT INE  EXCUTE ( ANUMB 1 , ANUMB 2, PO I  NT *  XA , XB , H I , NO , S T AR ) 
LOGICAL  SK IPlf SKIP2 
LOGICAL  LESS 

DIMENSION  XA (2000), XB *2000) 

INTEGER  POINT* 130, 100)*DASH,BLANK*STAR 

DATA  K1/10/*K2/100/*K3/ 1000/ *K4/l 0000/, NEG0NE/-1/ 

LESS-. FALSE. 

1-0 

J*0 

CALL  MEMORY*  ANUMB1*  ANUMB2, NO,  XA* XB I 
IF*  ABS (  ANUMB 1 ).GT. 10000.0)  GO  TO  50 
IF*  ABS*  ANUMB2 ) .GT. 10000.0)  GO  TO  50 
IF*  ABS*  ANUMB  1  I.EQ.O.O)  GO  TO  500 
IF* ABS* ANUMB21.EQ.0.0}  GO  TO  500 
IF*  ABS*  ANUMB2  ).LT. 0.001 )  GO  TO  50 
IF* ABS* ANUMB 1 )«LT« 0.001 )  GO  TO  50 

CALL  SCALE1(K1*K2*K3*K4*NEG0NE  *1 ,ANUM6l, ICONS, LESS, SKI  PI ) 
CALL  SCALE2(K1*K2*K3*K4* J, JCONS, ANUMB2* SKIP2 ) 

CALL  MPOINTtJ* JCONS* L*ANUMB2* SKI P2 ) 

CALL  SPOINT ( I  , ANUMBl* LESS, ICONS,  L* POINT , SKtPl,NO,xA,XB,STAR) 

RETURN 

END 


SUBROUTINE  SCALE1(K1,K2,K3,K4*NEG0NE, I, ANUMBl* ICONS, LESS'  SKI  PI) 
LOGICAL  LESS 
LOGICAL  SKIP1 
SKIP1-. FALSE. 

I  CONS* 1 

1-1 

AKEEP- ANUMBl 
NUMBi -ANUMB 1 
NUMB- 1 ABS* NUMB  1) 

IHNlHB.EQ.O'oOTO  12 
IF* NUMB1 131*40*40 
IFINUMBl.GE. lOIGOTO  11 
GO  TO  50 
LESS-.TRUE* 

XXX— NUMB  l 

75 


84  ] 
842 
B43 

844 

845 
646 
84  7 
840 
849 
b!>0 
851 
Hb? 
Hb  3 
854 
8bb 
8b6 

857 

858 

859 

860 
861 
862 

863 

864 

865 

866 

867 

868 

869 

870 

871 

872 

873 

874 
075 

876 

877 
H78 

879 

880 
881 
882 

883 

884 

885 

886 

887 

888 

889 

890 

891 

892 

893 

894 

895 

896 

897 

898 

899 

900 


I 


non 


NUHB1-XXX 
GOTO  40 

12  G0f0( l«2i  J«4)  •  1 

1  ANUMB1-AKEEP 
IC0NS-K1 
RKi-Xl 

ANUMB1-AMUMB1*RK1 

1*2, 

GOTO  21 

2  ANUMB1*AKEEP 
IC0NS-K2 
RK2-K2 

ANUMB1-ANUMB1*RK2 

1-3 

GOTO  21 

3  ANUMB1-AKEEP 
IC0NS-K3 
RK3-K3 

ANUMBI«ANUMB1*RK3 

1*4 

GOTO  21 

4  ANUMB1-AKEEP 
IC0NS-K4 
RK4-K4 

ANUMBX*ANUMB1*RK4 
NUHB1*ANUMB1 
IF(NUMBl.EQ.O)  GO  TO  51 
GO  TO  50 
11  SKIPi-.TRUE. 

GO  TO  (6«7»8f 9)» I 

6  ANUMB1-AKEEP 
IC0NS-X1 
RKl-Kl 

ANUMB1-ANUMB1/RK1 

1*2 

GOTO  21 

7  ANUMB1-AKEEP 
IC0NS-K2 
RK2-K2 

ANUMB1-ANUNB1/RK2 

1*3 

GOTO  21 

8  ANUMBl-AKEEP 
IC0NS-K3 
RK3-K3 

ANUN81-ANUMB1/RK3 

1-4 

GOTO  21 

9  ANUNB1-AREEP 
1C0NS-K4 
RK4-X4 

ANUHB1-ANUMB1/RK4 
I*-.  GO  TO  50 
r-  51  1-5 


901 

902 

903 

904 

905 

906 

907 

908 

909 

910 

911 

912 

913 

914 
9  15 

916 

917 
916 
9  19 

920 

921 

922 

923 

924 

925 

926 

927 

928 

929 
9  30 

931 

932 

933 

934 

935 

936 

937 

938 

939 

940 

941 

942 
94  3 

944 

945 

946 

947 

948 

949 

950 

951 

952 

953 

954 

955 

956 

957 

958 

959 

960 
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SUBROUTINE  SCALE2IK1* K2t K3»  K4« J» JCONS t ANUHB2  t SKI P2 ) 

961 

LOGICAL  SKIP2 

962 

SKIP2-. FALSE. 

96  3 

JCONS-l 

964 

J«l 

965 

ANUMB2*  ABS ( ANUHB 2 1 

966 

6KEEP-ANUHB2 

967 

20  NUMB2-ANUNB2 

968 

IF( NUMB2.EQ.0)G0  TO  10 

969 

IFCNUMB2.GE.10IG0  TO  11 

970 

GO  TO  50 

971 

10  GOTO! 1* 2*3.4) *  J 

9  7? 

1  ANUMB2-BKEEP 

9  7  J 

JC0NS-K1 

9  74 

RKl-Kl 

975 

ANUMB2-ANUMB2*RKi 

9  76 

J-2 

977 

GO  TO  20 

97H 

2  ANUMB2-BKEEP 

979 

JC0NS-K2 

9a0 

RK2*K2 

9H1 

ANUHB2- ANUMB2*RK  2 

982 

J*3 

983 

GO  TO  20 

984 

3  ANUMB2-BKEEP 

985 

JC0NS-K3 

986 

RK3-K3 

987 

ANUMB2-ANUMB24RK3 

9HH 

J-4 

989 

GO  TO  20 

990 

4  ANUMB2-BKEEP 

991 

JC0NS-K4 

992 

RK4-K4 

993 

ANUHB2- ANUHB2*RK4 

994 

NUMB2-ANUMB2 

995 

IFINUMB2.EQ.0)  GO  TO  51 

996 

GO  TO  50 

997 

11  SKIP2-.TRUE. 

998 

GO  TO  (6f7«8«9)t J 

999 

6  JCONS-Ki 

1000 

ANUMB2-BKEEP 

1001 

RKl-Kl 

1002 

ANUMB2- ANUMB2/RK 1 

1003 

J-2 

1004 

GO  TO  20 

1005 

7  JC0NS-K2 

1006 

ANUNB2-BKEEP 

1007 

RK2-K2 

1008 

ANUMB2-ANUHB2/RK2 

1009 

J-3 

1010 

GO  TO  20 

1011 

8  JC0NS-K3 

1012 

ANUHB2-BKEEP 

1013 

RK3-K3 

1014 

ANUMB2 ■ ANUtt  B2 /RK  3 

1015 

J-4 

1016 

GO  TO  20 

1017 

9  JC0NS-K4 

1018 

ANUMB2-BKEEP 

1019 

RK4-K4 

1020 

77 


ANUHB2*ANUHB2/RK4 

1021 

60  TO  50 

1022 

$1 

J«5 

1023 

50 

RETURN 

1024 

END 

1 0  2  *3 

C 

1026 

C 

1027 

C 

1028 

C 

1029 

SUBROUTINE  WPOINTI J* JC0NS*L,ANUHB2. SKIP2 » 

1030 

LOGICAL  SKIP2 

1031 

L-0 

1032 

I0ELTA-0 

1033 

IFUNUHB2.GE.9.2)  IDELTA«12 

10  34 

IF( ANUHB2.GE.0.O)  GO  TO  200 

1035 

IFUNUHB2.GE.6.9)  GO  TO  201 

10  36 

IFUNUHB2.GE.5.9)  GO  TO  202 

10  3  7 

IFUNUHB2.GE.5.0)  GO  TO  203 

10  3H 

IFUNUHB2.GE.4.2)  GO  TO  204 

10  39 

IFUNUHB2.GE.3.5)  GO  TO  205 

1040 

IFUNUHB2.GE. 2*9 )  GO  TO  206 

1041 

IFUNUHB2.GE.2.4 )  GO  TO  207 

1042 

IFUNUHB2.GE.1.9)  GO  TO  20B 

104  3 

IF(  ANUMB2.GE.U5)  GO  TO  209 

1044 

IFUNUHB2.GE.1.2)  GO  TO  210 

1045 

GO  TO  211 

1046 

200 

IFUNUHB2.LT.9.2 )  IDELTA-ll 

1047 

201 

I F f ANUMB2.LT. 8.0)  IDELTA-10 

1048 

202 

IF( ANUHB2.LT.6.9)  IOELTA-9 

1049 

203 

IFUNUHB2.LT. 5.9)  I0ELTA*8 

1050 

204 

I F 1 ANUHB2.LT. 5.0  )  IDELTA«7 

1051 

205 

I F 1 ANUHB2.LT. 4. 2  )  IDELTA-6 

1052 

206 

IF( ANUHB2.LT. 3.5 )  IDELTA»5 

1053 

207 

IF( ANUHB2.LT. 2.9 )  IDELTA-4 

1054 

208 

I F ( ANUHB2.LT. 2.4 )  IDELTA-3 

1055 

209 

IFUNUHB2.LT.  1.9)  IDELTA*2 

1056 

210 

1FI ANUHB2.LT. 1.5)  IOELTA-l 

1057 

211 

IFUNUHB2.LT.  1.2 )  I0ELTA-0 

1058 

I DELTA* 13-1 DELTA 

1059 

IF( J.EQ.5)  GO  TO  50 

1060 

IFISKIP2)  GO  TO  41 

1061 

IFI JCONS.EQ* 1 ) GOTO  1 

!.062 

IFI JC0NS.EQ.101G0T0  10 

1063 

IF( JCONS.EQ. 100) GOTO  100 

1064 

IFI JCONS.EQ. t 000) GOTO  1000 

1065 

IFI JCONS.EQ. 10000 )GOTO  IO000 

1066 

1 

L-73+I0ELTA 

1067 

GO  TO  40 

1068 

10 

L*86+I0ELTA 

1069 

GO  TO  40 

1070 

100 

L-994IDELTA 

1071 

GO  TO  40 

1072 

•too 

L*ll2+IDELTA 

1073 

r*- 

GO  TO  40 

1074 

fttooo 

L*125 

1075 

00 

GO  TO  40 

1076 

in  4l 

KC0NS*JC0NS+1 

1077 

SKIP2-. FALSE. 

107P 

IFIKCONS.EQ. 11 )GOTO  11 

1079 

1FIKC0NS.EQ. 101 ) GOTO  101 

1080 

78 


58780 


Q  IFIKCONS.EQ.lOOl )GOTO  1001  1081 

IFI KCONS»EQ« 1000 1 ) GOTO  10001  1082 

11  L«60+I0ELTA  1084 

GO  TO  40  1084 

101  L-47MDELTA  1085 

GO  TO  40  1086 

1001  L»34*IDELTA  1087 

GO  TO  40  1068 

50  L«125  1089 

GO  TO  40  1090 

10001  L*34  1091 

40  RETURN  1092 

END  1094 

C  1094 

C  1095 

C  1096 

C  1097 

SUBROUTINE  SP01NTI ! , ANUMB l, LE S S, ICONS, L, P01 NT , SKI  PI ,N0 , X A fXH , STAR )  1096 

DIMENSION  XA( 2000 ),XB(2000)  1099 

INTEGER  P0INT(130,100),DASH,BLANK,STAR,Q  1100 

INTEGER  STARl,STAR2  1101 

DATA  STAR1/1H0/,STAR2/IH*/  1102 

DATA  BLANK/ 1H  /,0ASH/1H-/,  Q/lHI/  1103 

LOGICAL  SKI  PI  (LESS  1104 

AKEEP»ANUMB1  1106 

ANUMBl* ABSI AKEEP I  1106 

K00L*0  1107 

IFU.EQ.5IG0  TO  1000  1106 

IFIANUMB1.GE.7.1IK00L-7  1109 

IFIANUMB1.GE.5.0IG0  TO  1011  1110 

IFI ANUMB1.GE.4.01GO  TO  1111  1111 

I FI ANUMBI.GE.3. 1 1 GO  TO  1211  111? 

IFI ANUMBI.GE.2.3IG0  TO  1311  1113 

IFC ANUMB1.GE.1.6IG0  TO  1411  1114 

IFIANUMBl.GE.l.OIGO  TO  1511  1115 

GO  TO  4011  1116 

1011  IFI  ANUMBl.LT.7.1  IKOOL-6  1U7 

1111  IFIANUMB1.LT.5.0)K00L*5  1118 

1211  IFI ANUMB1.LT.4.0IK00L-4  1119 

1311  IFI ANUMB1.L  T.  3. 1 1K00L *3  1120 

1411  IFI ANUMB1.LT.2.3IK00L-2  1121 

1511  IFI ANUMB l.LT. 1*6 )K00L*1  1122 

4011  NCRMNT-8-K00L  1123 

IFILESSIGO  TO  40  1124 

IFISKIPDGO  TO  41  1125 

IFI IC0NS.EQ.1IG0T0  1  1126 

IFI ICONS. EQ. 10IG0T0  10  1127 

IFI ICONS.EQ.IOOIGOTO  100  1128 

IFI ICONS.EQ.IOOOIGOTO  1000  1129 

1  KOOL-7WKOOL  1130 

GO  TO  50  1131 

10  K00L-644K00L  1132 

GO  TO  50  1133 

100  KOOL-574KOOL  1134 

GO  TO  50  1135 

000  KOOL-504KOOL  1136 

GO  TO  50  1137 

41  LCONS-ICONS+l  1138 

SKIPl-.FALSE.  1139 

IFILC0NS.EQ.11IG0T0  11  114,0 

79 


I 


5878J 


c 

B 


1 F! LCONS*  EQ« 101 1G0T0  101  1 1  «♦  1 

IF  I LCONS*  EQ« 1001 (GOTO  1001  l  LA? 

11  K00L=78*K00L  11A3 

GO  TO  50  l  144 

101  K00L-85+K00L  1146 

GO  TO  50  1146 

1001  K00L«924K00L  1147 

50  LESSs«FALSE,  H46 

IF!POINT!L,KOOL).EQ.DASH)  GO  TO  51  1147 

IF (  POINT ( L»KOOL ) «EQ«0  1  GO  TO  51  1156 

14  IFIPOINTlL.KOOU.EQ.STARl)  GO  TO  15  1151 

IF!P0INT!L,K00U.E0.STAR2)  GO  TO  15  115? 

IF! STAR* EO* STAR  1 )  GO  TO  51  1163 

IFISTAR.EQ.STAR2)  GO  TO  51  1 1 

IF!  POINT! L«KOOL )«NE* BLANK )  GO  TO  70  1155 

51  POINT ( L«KOOL 1 *STAR  1156 

ANUMB l* AKEEP  1157 

GO  TO  70  .11 5a 

15  L*L  ♦  l  1159 

IF(L.E0*130)  GO  TO  51  1160 

GO  TO  14  llol 

40  IF! SK (PI  )  GO  TO  42  lib? 

IZ«IC0NS42  llbi 

IF! IZ.EQ, 31G0T0  3  1164 

IF!  IZ • E0« 1 2 ) GOTO  12  L  165 

IF!  IZ.E0.1021G0T0  102  1166 

IF!  IZ.EQ.1002IG0T0  1002  1167 

3  K00L*2 l ♦NCRMNT  1165 

GO  TO  50  1164 

12  K00L*28+NCRMNT  11 70 

GO  TO  50  U  71 

102  K00L-354NCRKNT  117? 

GO  TO  50  117  J 

1002  KOOL* 42 ♦NCRMNT  1174 

GO  TO  50  1175 

42  IYMCONSO  1176 

SKIP1*.FALSE.  U77 

IF! IY.EO. 131G0T0  13  1 178 

IF!  IY.EO. 1031G0T0  103  1174 

IF!  IY.EO.  10031G0T0  1003  UuO 

lF( IY.EO. 10003 1G0T0  10003  lldl 

13  K00L»144NCRMNT  11k? 

GO  TO  50  1 1 n  3 

103  K00L«7^!CRMNT  1164 

GO  TO  50  1185 

1003  KOOL-NCRMNT  lltj6 

GO  TO  50  1187 

10003  POINT ! L , l ) *STAR  1168 

GO  TO  70  1189 

70  RETURN  1190 

END  1191 

C  119? 

1193 

1194 

1195 

SUBROUTINE  PREPARIPOINT »  1196 

INTEGER  POINT!  130, 1001, PRT  1197 

OATA  PRT  /  b/  1198 

WR I TE ! 6, 1 1  1149 

1  FORMAT! 1HI  1  1200 


80 


WRITE!6,11) 

1201 

11  FORMAT!///) 

1202 

WRITE (PUT* 14) 

1203 

14  FORMAT !75X,8HL0G  PL0T,/58X,42HC0MPLEX  FREQUENCY 

PLANE, LEFT  HANC  CU 

1204 

*ADRANT,/7SX,9H!RAD/SEC)  ) 

1205 

WRITE16, 12) 

1206 

12  FORMAT!  32X,5H10000,9X,4HIOOO,9X,3H1CO,10X,2H10, 

12X, IH1 , 1 IX, 2H, 1,10 

1207 

•X*  14H,01< - J-OMEGA, /  18X,11HMINUS  SIGMA) 

1208 

13  FORMAT! 32X,5HlOOOO,9X,4HlOOO,9X,3HiOO,lOX,2HlO, 

12X, 1H1, 1 IX , 2H* 1,  10 

1209 

*X,14H.0l< - J-OMEGA*/  IMX.11H  PLUS  SIGMA) 

1210 

00  SO  1-1,130 

1211 

50  P0INT!I,50)«P0lNT!l,r>l) 

1212 

WRITEIPRT, 10) POINT 

1213 

10  FORMAT! IX, 130A1) 

1214 

WRITECo, 13) 

1215 

WRITEIPRT, 15) 

1216 

15  FORMAT! 75X,8HL0G  PLOT,/58X,43HCOMPLEX  FREQUENCY 

PLANE, RIGHT  HANO  0 

1217 

*UADRANT,/75X,9H| RAO/SEC)  ) 

1 2 1H 

RETURN 

1219 

ENO 

1220 

c 

1221 

c 

1222 

c 

1223 

c 

1224 

SUBROUT I NE  MEMOR Y ! ANUMB 1 • ANUMB 2,  NO , X  A  ,  X  B ) 

1225 

DIMENSION  XAI2000),XB!2000) 

1226 

XAINO l-ANUMBl 

1227 

XB! NO l-ABS! ANUMB 2 ) 

1228 

NO«NO*l 

1229 

RETURN 

1230 

ENO 

1231 

-  c 

1232 

i  C 

1233 

C 

1234 

C 

1235 

SUBROUTINE  WRITI T!XA, XB ) 

1236 

DIMENSION  XA!2000),XB!2000) 

1237 

WR ITE! 6, l ) 

1238 

1  FORMAT! 1H1) 

1239 

WRITE!*, Ill 

1240 

11  FORMAT! 5X,48HTHE  FOLLOWING  ROOTS  ARE  PLOTTEO  ON  THE  LOG  PLOT,/ 

1241 

15X,99HR00TS  AT  THE  ORIGIN  ARE  NOT  PRINTED  OR  PLOTTEO,  ROOTS  ON  THE 

1242 

•  J-OMEGA  AXIS  ARE  NOT  PLOTTEO, 

, 

1243 

2, //16X,5HSIGMA,?SX,7H J-OMEGA,//) 

1244 

00  77  IZO-1,2000 

1245 

XC-XB1 120) 

1246 

IF!XC,GT,0,0)  XC— XC 

1247 

IF! XC,NE,0,0)  GO  TO  22 

1248 

20  IF! XAI IZOI 122,40,22 

1249 

40  XI- IZO 

1250 

K2-UO+12 

1251 

00  50  K-K1,K2 

12S2 

CVJ 

IF! XAIK ),NE,0,0)  GO  TO  77 

12S3 

00 

SO  IF!XB!K),NE,0,0)  GO  TO  77 

1254 

r- 

GO  TO  60 

1255 

00 

22  WRITE!6,10IXA!IZ0),XC 

1256 

10  FORMAT! SX,F20,9, 10X,F20S9) 

1257 

• 

77  CONTINUE 

1258 

60  RETURN 

1259 

♦—  » 

END 

1260 

« 

81 

Security  Clsssificstion 


DOCUMENT  CONTROL  DATA  •  R  &  D 

Security  clmtai  fie  atton  of  tltlm,  body  oi  abatract  mnd  indexing  annotation  muat  bo  antarad  whan  tha  overall  raport  la  etaaalllad) 


l  o-i  NG  A  C  T I V  l  T  V  (Corporal*  author)  1 2a.  REPORT  IKCURI  T  Y  CLASSIFICATION 

U  S.  Army  Aberdeen  Research  and  Development  Center  |  Unclassified 
Army  Materiel  Systems  Analysis  Agency 
Aberdeen  Proving  Ground,  Maryland 


3  REPORT  TITLE 

\  i  inear  Closed  Loop  System  Analysis  Procedure  Using  Line  Printer  Plots  of 
Characteristic  Equation  Root  Loci 


|a».  group 


4  DESCRIP  TIVE  NOTES  (Typa  of  rape* I  and  Include  dalaa) 


t  * u  T mi>  R IS)  (Flrat  nama,  ailddta  Initial,  fast  nama) 

Harold  H.  Burke,  Robert  L.  Payne,  Jr. 


fl  MCPOR  T  DATI 

November  19o8 


i.  CONTRACT  OR  GRANT  NO. 


7a.  total  no  of  RAGES 

85 


M.  ORIGINATOR'S  REPORT  NUMStRISI 


1b.  NO.  OF  REFS 


»>  project  no.  RDTt,E  lT523801A0'Jb 


Technical  Memorandum  No.  21 


s b.  OTHER  REPORT  noisi  (Any  otnar  numbara  dial  mtay  ba  aaalgnad 
tfils  raparl) 


io  oistrigution  statement 


Tins  document  is  subject  to  special  export  controls  and  each  transmittal  to  foreign 
governments  or  foreign  nationals  may  be  made  only  with  prior  approval  of  Commanding 
Officer.  U.  S.  Army  Aberdeen  Research  and  Development  Center,  Aberdeen  Proving  Ground 

M  I  V»*  I  RI»<1  “ 


4  M  "VO  MFtACIIOOROffy  I4TI, 
VGA*"#  /  O  OHOkin  FOR  ARMY  USB 


OO  FORM  I4TS.  •  JmM  M.  WHICH  II 


Unclassi fied 

Swcurily  Clsssihcatiofl 


Root  Locus 
Plot  Routine 

linear  Control  Stability  System  Analys 


